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Abstract. Ellerman bombs are bright emission features observed in the wings of 
Ha, usually in the vicinity of magnetic concentrations. Here we show that they can 
also be detected in the Ca II infrared triplet lines, which are easier to interpret and 
therefore allow for more detailed diagnostics. We present full Stokes observations of 
the 849.8 and 854.2 nm lines acquired with the new spectro-polarimeter SPINOR. 
The data shows no significant linear polarization at the level of 3x 10 -4 . The circular 
polarization profiles exhibit measureable signals with a very intricate pattern of 
peaks. A non-LTE analysis of the spectral profiles emerging from these features 
reveal the presence of strong downflows (~10 km/s) in a hot layer between the 
upper photosphere and the lower chromosphere. 
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1. Introduction 

One of the most intriguing discoveries made in Solar Physics early last 
century was the phenomenon known as Ellerman bombs. The common 
use of the Ha line to monitor the Sun allowed the discovery of features 
with enhanced emission peaks in the wings of this line, not affecting 
the central absorption core (Ellerman 1917). Off-band (~ 1 A) Ha 
filtergrams showed these features to be spatially point-like structures 
with a size of lasting for little more than 10 minutes (see Qiu 

et al. 2000 and references therein). As they are virtually absent from 
line core filtergrams, they have always been associated to events taking 
place in the low chromosphere. This is also confirmed by their clear 
correlation with brightenings seen in the 1600 A continuum images from 
the TRACE satellite (Qiu et al. 2000). Sometimes, an X-ray transient 
brightening can also be observed at the same location (Shimizu et al. 
2002). Recent reviews on Ellerman bombs are given by Rust (2001) 
and in the introduction of Georgoulis et al. (2002). 
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The physical process behind Ellerman bombs has remained elusive. 
The blue asymmetry observed by the Ha profiles (excess emission in 
the blue wing) has sometimes been associated to outward motions of 
the solar material (Koval and Severny 1970). But using two-component 
models separated along the line of sight, Dara et al. (1997) showed that 
the asymmetry could be originated from, either, hot emitting material 
moving upwards or cold absorbing material falling from the top of 
the chromosphere. The relation to the magnetic field configuration has 
been controversial as well. Ellerman bombs have been detected in two 
distinct solar magnetic scenarios, but always inside active regions. They 
occur in emerging regions between the opposite polarity footpoints and 
also surrounding isolated sunspots, beyond the penumbral boundary. 
For this later case, Nindos and Zirin (1998) have shown that they can 
be identified with the so-called Moving Magnetic Features (MMFs) 
that are seen in the moat region around sunspots (for a description of 
MMFs see Harvey and Harvey 1973). Whereas this points towards a 
magnetic origin of Ellerman bombs, a number of events observed by 
Nindos and Zirin (1998) (see also Dara et al. 1997) were not clearly 
linked to magnetic structures. These authors point out that a number 
of cases showed a relationship with the boundaries of magnetic regions 
but others did not. 

More recently the situation has been clarified by the analysis of the 
observations made by the Flare Genesis balloon experiment (Georgoulis 
et al. 2002, Pariat et al. 2004). Observing an emerging flux region, 
these authors describe dipolar features that appear in the middle of 
the region, where the flux emerges at the surface (below the so-called 
arch filament system). With better sensitivity magnetograms than the 
previous works, Georgoulis et al. (2002) found that almost all Ellerman 
bombs studied could be identified with the neutral line above the dipo- 
lar features (or in neutral lines inside supergranular boundaries). Those 
cases not linked with the presence of neutral lines could be associated 
with interacting field lines above the photosphere, a result confirmed 
by the extrapolations made by Pariat et al. (2004). These authors 
proposed, then, as the most probably mechanism for the generation 
of Ellermam bombs some form of magnetic reconnection in the low 
chromosphere. That reconnection happens efficiently at these heights 
is supported by the fact that the temperature minimum is the region 
with the lowest electrical conductivity, and resistive processes occur 
there more easily than elsewhere (Litvinenko 1999). While the Flare 
Genesis observations correspond to the central portion of an emerging 
active region, these conclusions offer an interesting parallelism with 
the Ellerman bombs observed surrounding sunspots and associated to 
MMFs. Pariat et al. (2004) pointed out that many of the bombs studied 
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were spatially located at the dipped portions of undulatory field lines 
in the flux emerging area. This undulatory or serpentine configuration 
is also used for the explanation of MMFs (as first proposed by Harvey 
and Harvey 1973). In this way, a unifying scenario for the generation 
of these point-like energy release processes could be constructed. But 
a study similar to that made by Pariat et al. (2004) for the moat of 
sunspots is still lacking. The comparison between the Ellerman bombs 
seen in these two scenarios and the determination of their physical 
origin awaits for more solid observational tools than Ha line profiles. 

In particular, techniques providing information of the magnetic field 
at the chromospheric heights where the energy release takes place, will 
be of great value. Ha linear polarization signals observed in Ellerman 
bombs have been reported in the past (Rust and Keil 1992), but their 
exact amount an origin remains controversial. In this paper, we propose 
the use of full Stokes spectropolarimetry in more reliable diagnostic 
lines such as the Call infrared triplet lines. The most important advan- 
tage these lines offer is the availability of powerful inversion techniques 
that can be used to cover almost the full chromospheric region. This 
paper presents a first characterization of the atmosphere of two regions 
close to a sunspot that display the properties ascribed to the Ellerman 
bombs. 



2. Observations 

The observations used in this work are first-light data from the new in- 
strument Spectro-Polarimeter for INfrared and Optical Regions (SPINOR, 
see Socas-Navarro et al. 2005), at the Dunn Solar Telescope (DST) 
located at the Sacramento Peak observatory. The dataset contains the 
two chromospheric Call lines at 849.8 and 854.2 nm and two photo- 
spheric Fel lines at 849.7 and 853.8 nm. We observed active region 
NOAA 0634 on June 16 at 15:16 UT, carrying out a 350-step scan with 
a spatial stepping of 0.22". By using the adaptive optics system of the 
DST (Ren et al., 2003), we attained a spatial resolution as good as 0.6" , 
although this figure varies slightly along the scanning direction (the 
horizontal direction in the figures) due to temporal seeing variations. 

In order to improve the signal-to-noise ratio in the Stokes profiles, we 
averaged the observations in 3x3 pixel boxes. The noise is thus reduced 
to ~3xl0 -4 times the continuum intensity. At this level there is barely 
no signal in Stokes Q and U, so we decided to leave them out of our 
analysis. The chromospheric Stokes V profiles, however, reveal fasci- 
natingly complex shapes with multiple peaks, while the photospheric 
lines exhibit the usual antisymmetric profiles. 
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Figure 1. Top: Continuum map reconstructed from the spectro-polarimetric obser- 
vations. North is up and west is left of the image. The arrow indicates the direction 
to the solar limb. Bottom: Circular polarization in the red lobe of the photospheric 
Fe I line at 8497 A. The horizontal dark stripes are hairlines inserted in the slit to 
serve as a reference for the calibration. 
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Figure 2. BBSO Ha image of NOAA 0634 taken on 16 June 2004 at 15:50 UT. 
Notice the bright emission on the western side of the sunspot, near the penumbra. 



Fig 1 shows maps of the continuum intensity and a "magnetogram" 
taken in one of the photospheric lines. The magnetic configuration 
surrounding the sunspots appears very different on both sides of it. The 
eastern side exhibits a much more complex topology with an intricate 
pattern of plage fields extending from the sunspot boundary to the edge 
of the field of view. Ha images of the region (Fig 2) reveal considerable 
emission near the eastern penumbral boundary. 
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2.1. Velocities 



Photospheric velocity maps of NOAA 0634 were obtained by measuring 
the minimum position of the Fe I line at 849.7 nm. A similar strategy 
is not viable for the Ca II lines, however, because of the complicated 
patterns of emission reversals and self-absortions found in the line cores. 
Instead, we measured the intensity difference between two points sym- 
metrically located at various distances from the line core. This method 
is only an approximation but it works better than taking the minimum 
position. 

Fig 3 shows the velocity maps obtained at two heights, with the 
convention that positive velocities are away from the observer. The 
cosine of the heliocentric angle for this region is jjl = 0.71. Notice the 
two prominent bright features marked with arrows in the lower panel 
of the figure. These features are strong plasma flows directed away 
from the line of sight and located on the center side of the sunspot 
(but outside the photospheric penumbral boundary). Detailed profiles 
of these features are shown in Fig 4. 

Both features exhibit very strong redshifts, visible in Stokes I and V. 
Point A shows a very prominent redshift starting near 40 pm away 
from line center (corresponding roughly to the high photosphere or 
low chromosphere) all the way to the line core (middle chromosphere) . 
This feature is well localized (see lower panel) and the Doppler shift is 
~14 km s" 1 . Point B exhibits somewhat weaker flows (~10 km s _1 ) that 
are more localized in height (upper photosphere/lower chromosphere), 
but has a more extended tail along the slit (see lower panel). Notice 
how both points A and B coincide with Ha emissions in Fig 2 



3. Non-LTE inversions 

We performed detailed non-LTE inversions of the area surrounding 
points A and B in Fig 3. The inversions were carried out with the code 
of Socas-Navarro et al. (2000b) (see also Socas-Navarro et al. 1998) 
by applying it systematically on each pixel independently. In order to 
minimize the risk of the algorithm settling on a secondary minimum, 
we repeated each inversion 10 times with randomized initializations. 
The model atmosphere that yields the best fit to the observed Stokes 
profiles is taken as a good representation of the solar atmosphere under 
study. 

The inversion code treats the Call lines in non-LTE and the Fel 
lines in LTE. Blends are treated consistently by adding the opacities of 
all the relevant transitions at each wavelength and spatial point, and 
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Figure 3. Top: Velocity map (arbitrary units) measured as the shift of Fe I AA8497 
line core. In this paper we follow the Astrophysical convention that positive values 
correspond to redshifts. Bottom: Dopplergram obtained by subtracting the intensity 
at two symmetric positions of the Ca II line at 8498 A. The positions measured are 
1.35 A apart, on the wings of the Ca II line. The photospheric character of both 
maps is evident from the granulation pattern. The areas marked by arrows and 
labeled "A" and "B" in the figure exhibit strong (~10 km s _1 ) flows and have been 
subject to detailed analysis in this work. 
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Figure 4- Top: Two sample profiles of the Ca II line at 8498 A observed in areas 
A (left) and B (right). These profiles exhibit particularly striking downflows. The 
dashed line shows the fit from the non-LTE inversion. Overplotted in dotted line for 
reference is the average profile along the slit direction. Bottom: Stokes V spectra 
along the slit. The arrows mark the position of the profiles shown in the upper 
panels. 



then solving the radiative transfer equation with the total opacity. The 
non-LTE calculation is based on the preconditioning strategy of Ry- 
bicki and Hummer (1992) (see also Socas-Navarro and Trujillo Bueno 
1997) assuming complete angle and frequency redistribution. This last 
approximation is particularly good for the Call infrared triplet lines, 
as shown by Uitenbroek (1989). 

Let us consider in detail the models obtained for the points that we 
labeled as A and B. Figures 5 and 6 show the atmospheric stratification 
inside the downflows, as well as an average over the surrounding, non- 
downflowing area for reference. The most striking features in these 
models are the strong downwards velocities found in the upper pho- 
tosphere, around log(rsoo) = —3. The motions along the line-of-sight 
approach the speed of sound in point B and exceed it in point A. 
The other parameters show different behaviors in points A and B. In 
particular, the downflowing layer in point B is considerably hotter than 
the reference atmosphere. However, in point A the downflowing layer is 
cooler (although it then becomes hotter at higher layers). In both cases 
there is another hot layer at log(Tsoo) = —5. Our findings for the two 
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Figure 5. Average model atmosphere in the downflowing feature labeled as point A 
(solid) and its surroundings (dashed). Panels from top to bottom, left to right 
are temperature, line-of-sight velocity, magnetic field strength and inclination. The 
dotted line in the upper right panel represents the speed of sound in the atmosphere. 
Error bars are plotted at the inversion nodes (see Socas-Navarro et al 2000; 2001 
for details). 



regions are consistent with the results of Georgoulis et al. (2002) who 
found Ellermam bombs to be predominantly dominted by downflows 
with no indication of upflows. 

The fits provided by the inversions can be seen in Fig 7. Note that the 
observed chromospheric profiles exhibit a very complex structure with 
multiple peaks in Stokes V (the photospheric lines, on the other hand, 
produce rather normal profiles). Generally speaking, the inversion code 
yields a reasonably good fit to the overall features, but there is more 
structure in the observed profiles (particularly noticeable in Stokes V 
with multiple peaks), than in the synthetic ones. This is an indication 
that the actual atmospheres are even more complex than our models, 
probably due to the presence of unresolved atmospheric components. 

Thus far, we have not made a direct connection between the down- 
flowing features analyzed here and EBs. To this aim, we used the model 
atmospheres obtained from the inversion to synthesize the Ha profiles 
that these models produce. The result is plotted in Fig 8, along with a 
reference quiet Sun profile calculated from model C of Vernazza et al. 
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Figure 6. Average model atmosphere in the downflowing feature labeled as point B 
(solid) and its surroundings (dashed). Panels from top to bottom, left to right 
are temperature, line-of-sight velocity, magnetic field strength and inclination. The 
dotted line in the upper right panel represents the speed of sound in the atmosphere. 

(1981) (VAL-C). This synthetic profile has the same appearance tradi- 
tionally observed in EBs. Ellerman (1917) describes the Ha spectra of 
his bombs as having emission features (brighter than the continuum) 
on both sides of the line center. The red asymmetry (higher intensity 
peak in the red wing) shown here is not the one most commonly ob- 
served in Ellerman bombs (see the introduction), but this asymmetry 
is also present quite often (40% of the cases according to Rust 2001). 
It is the result of a hot downflowing component as seen in Fig 5. The 
resemblance of EBs with the synthetic profile in Fig 8 is a solid evidence 
that we are indeed observing the same phenomenon. Moreover, it also 
provides strong support for the reliability of the inversions. 

4. Conclusions 

The energetic phenomenon generally referred to as EBs has remained 
very ellusive to solar physicists for almost a century. This is probably 
due in no small measure to the fact that it has only been observed in 
Ha thus far, which is a very difficult line to interpret. In this paper we 
present the first observations of EBs in other lines, namely the Call 
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Figure 7. Fits to 849.8 and 854.2 nm (left and right side of each panel, respectively). 
Left and right panels represent typical profiles in the downflowing areas in the vicin- 
ity of points A and B respectively. Upper panels: Stokes I. Lower panels: Stokes V. 
In all cases the solid (dashed) line represents the observations (fits). The profiles in 
this figure have been averaged over a 3x3 pixel box to improve the signal-to-noise 
ratio in Stokes V. 



infrared triplet. These lines are much easier to model and have been 
successfully inverted in the past few years (e.g., Socas-Navarro et al. 
2000a; Socas-Navarro et al. 2001; Socas-Navarro 2005b) and thus our 
results open a new window to the observational study of EBs. 

We find that the EBs analyzed in this work are associated with su- 
personic (or nearly so) downflows in the upper photosphere. No upflows 
are observed in our data. Socas-Navarro (2005a) determined the vector 
current densities in the large sunspot of Fig 1 and found that the east 
penumbra is permeated by strong currents in the upper photosphere, 
with a magnitude of ~3xl0 5 A km~ 2 . The fact that these EBs occur 
in a region with strong currents suggests that EBs may be associated 
with magnetic reconnection. 

In order to understand the nature of EBs and the physical processes 
that produce them, it would be extremely valuable to have time-series 
observations of photospheric and chromospheric lines, ideally with full 
Stokes polarimetry. 
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Figure 8. Synthetic Ha profiles obtained with the model atmosphere from the in- 
version of point A (solid) and the VAL-C model (dashed). Point B produced a 
qualitatively similar profile. 
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